Phase fluctuations in superconductors: 
from Galilean invariant to quantum XY models 
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We analyze the corrections to the superfluid density due to phase fluctuations within both a 
continuum and a lattice model for s- and d-wave superconductors. We expand the phase-only action 
beyond the Gaussian level and compare our results with the quantum XY model both in the quantum 
and in the classical regime. We find new dynamic anharmonic vertices, absent in the quantum XY 
model, which are responsible for the vanishing of the correction to the superfluid density at zero 
temperature in a continuum (Galilean invariant) model. Moreover the phase-fluctuation effects are 
reduced with respect to the XY model by a factor at least of order l/(fcF^o)'^. 



The issue of explaining the Unear temperature depen- 
dence of the superfluid density ps in cuprate supercon- 
ductors renewed the interest on the effects of the phase 
fluctuations (PF) of the order parameter in the supercon- 
ducting phase at low temperature, fndeed, both quasi- 
particles j^] and PF [|j in the classical limit have been 
suggested as alternative explanations for the observed de- 
pletion of Ps- Recently Refs. |^,|| considered the quan- 
tum effects on PF and showed that the classical limit is 
reached at a temperature Td which is too high to account 
for the low-temperature linearity of ps in cuprates. 

A crucial point in the analyses of PF effects is the 
choice of the effective model for the PF. In Refs. the 
Gaussian action for the phase was derived microscop- 
ically, while the anharmonic (non-Gaussian) terms were 
obtained by expanding the cos(6',; — 6j) coupling term of 
a lattice quantum XY model, derived by coarse-graining 
the Gaussian action on the scale of the coherence length 
^0, in powers of {6i — Oj) ^ 'Co|V0|. In Ref. Q, instead, 
also the anharmonic terms were derived microscopically, 
within a d-wave continuum BCS model. In this approach 
the interaction vertices for \I0 are determined by the the 
fermion loops and are smaller than the corresponding ver- 
tices of the XY model, leading to one- loop corrections to 
Ps which are much smaller than those of the quantum XY 
model, both at T = and at T > 0. However, they found 
a finite correction to ps even at T = 0, which is expected 
within a lattice XY model, but not within a continuum 
(Galilean invariant) model where ps equals the particle 
density p at T = This indicates that something is 
missing in their analysis. Moreover their description of 
PF effects within a continuum model is too restrictive, 
as it neglects lattice effects, which are certainly relevant 
in cuprates. 

In this paper we give a detailed analysis of the one- 
loop correction to ps due to PF within both a continuum 
and a lattice model for s- or c?-wave pairing. We specifi- 
cally consider the weak- to intermediate-coupling regime, 
even though at the end we also comment on the strong- 
coupling limit. We find that the microscopic derivation 
of the phase-only action, besides the classical (static) an- 
harmonic terms {\/6Y considered in Ref. p], introduces 



new third- and fourth-order quantum (dynamic) interac- 
tion terms which contain the time derivative of 9. These 
quantum terms are absent in the quantum XY model, 
where the dynamics only appears at the Gaussian level, 
and induce a correction to ps which cancels exactly, in 
the continuum case, the contribution due to the classical 
interaction, restoring the equality ps — p&tT^Q. On 
the other hand, the same cancellation does not hold in 
the lattice case, in which ps equals the average kinetic 
energy at T = 0. The inclusion of both classical and 
quantum interaction terms leads to a finite one-loop cor- 
rection to Ps, which is however of order l/(fci?^o)^ with 
respect to the result of the quantum XY model, kp be- 
ing the Fermi wave vector. In the classical regime we 
find that the PF correction to ps is smaller than within 
the classical XY model by the factor ^ I/(fci?^o)^> for 
both the continuum and lattice model. The reduction 
of the PF effects for fcj^^o ^ 1 is made even more pro- 
nounced in the continuum case by the inclusion of long- 
range Coulomb forces. The fact that the XY model 
generically overestimates the PF effects further supports 
the claim of Refs. Hl^, in which, even adopting the XY 
model, it has been shown that the contribution of PF 
does not account for the linear temperature decrease of 
Ps in cuprates and gives temperature-dependent small 
corrections when compared to the experimental data. 
We start with a continuum BCS action at a tempera- 



ture T — [3 ^ in d dimensions, S — J^^ drC, with 
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Hi - -nEk.k',q7k7k'C+_^a|C+k+a|C_k'+axCk'+af 
Here VI is the volume, > is the pairing interac- 
tion strength and the factor 7k controls the symmetry 
of the order parameter. In the following, unless explic- 
itly indicated, we set h = kg — 1. We perform the 
standard Hubbard-Stratonovich decoupling of Hj and 
make the dependence on the phase 6 of the complex 
order parameter A = |A|e*^ explicit by means of the 
gauge transformation Ca- Co-e*^/^ Then, af- 

ter integrating out the fermions around the supercon- 
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ducting saddle-point solution, and neglecting the fluc- 
tuations of |A| we obtain the effective action for 
PF, Seff[e] = TrJ2n=i ^(SGo)". Go is the mean-field 
Nambu Green function, and the self-energy matrix is 
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(V^- V)to, 
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where are the Pauli matrices, and the operator V= (V 

— V)/2 acts on Go- In the following we distinguish the 
"bosonic" contributions, generated by the T3 term in (|^), 
analogous to those of a boson model in the presence of 
condensate, the "fermionic" contributions, generated by 
the To term, and the "mixed" ones, obtained by combi- 
nations of To and T3 terms. Both bosonic and fermionic 
terms contribute to the Gaussian phase action which, in 
the hydrodynamic limit, reduces to the well-known form 



(3) 



where q = (q, x is the compressibility and D{T) = 
fi^Ps/m is the superfluid stiffness (in dimensional units). 
At Gaussian level the temperature dependence of D is en- 
tirely due to the quasiparticle excitations, giving D{T) = 
™ + l^Tv^^fiE^)- Here f'{x) = -/3e^V(e''- + 1)^, 
-^k = i/^k + is the quasiparticle dispersion in the su- 
perconducting state, = k^/2m— ^ is the free-electron 
dispersion, and Ak = |A|7k is the superconducting gap. 
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The one- loop PF corrections to D{T = 0) are induced 
by the anharmonic terms in 9. As we discussed above, 
within the quantum XY model the dynamics is intro- 
duced only at mean-field level adopting the PF Gaussian 
propagator P{q) = A{xtol + Dq^)-^ (see Eq. (|) |,|, 
and the one- loop correction to D, coming from the purely 
static interaction term cx S^qD J2a=ii9aS)^ , is 



SDxY = - 



Deo T 

2 dn 
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where g = C''+^/(i((i -f- 1), and C, ~ l/^o is the PF mo- 
mentum cutoff [0 . 

By contrast, we derive here the anharmonic terms by 
expanding Seff up to fourth order in 6 and evaluate 
the one-loop corrections to the Gaussian action. All the 
fermionic and mixed terms generate one-loop corrections 
to D which vanish or cancel each other at T = 0, what- 
ever is the symmetry of the gap, in agreement with the 
result previously discussed in Ref . Q , where however the 
contribution of the bosonic diagrams was incorrectly eval- 
uated. Indeed, the bosonic 6 term in Eq. (H) introduces 
a dynamic contribution in Sef / also beyond the Gaussian 
level, as it is represented by the third- and fourth-order 
bosonic vertices depicted in Fig. la. In particular, the 
presence of a third-order quantum interaction vertex, ne- 
glected in Ref. Q , leads to separate cancellations between 
both fermionic or mixed and bosonic contributions. Let 
us consider the one-loop diagrams for the bosonic cor- 
rections in Fig. lb: the first two are the self-energy 
corrections to Go, which shift the particle density p in the 
presence of PF, at fixed fj,. Alternatively, at fixed den- 
sity, we take these contributions into account by shifting 
/i with respect to its mean-field value, in order to keep 
p fixed. As a consequence, the only corrections to D 
with respect to p/m due to PF come from the last three 
diagrams in Fig. lb, i.e. 
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FIG. 1. (a): Bosonic diagrams in the non-Gaussian 
phase-only action up to fourth order. The dashed line in- 
dicates a 6 insertion, the solid hne represents Go. The vertex 
with a single incoming line corresponds to the insertion of 
the (6/2)ts term of Eq. (^, the vertex with two incoming 
lines to the insertion of ((V0)^/8m)r3. (b): Bosonic one-loop 
corrections to the PF Gaussian action. 

In the present context, to make comparison with Ref. 
1^ , we do not consider the effect of a dissipative term in 
Eq. (j^) for a c?-wave superconductor, even though this 
term would be essential for the low-temperature behavior 
and for the estimate of the temperature Td [|| . We first 
consider a neutral system, while the long-range Coulomb 
forces will be considered later. 



where xil) is the density-density bubble, which gives the 
compressibility x in the limit = 0, q — > 0. Notice that 
in writing Eq. (^) we are relying on the fact that both the 
third- and fourth-order vertices needed to calculate SD 
are expressed in terms of xil)- This is the relevant con- 
sequence of the Galilean-invariant form of the bosonic 
term appearing in the self-energy (^) |^Jl^. Evaluating 
Eq. (^ in the hydrodynamic limit, we get 
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dm^fl 
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where Eq = y^D/x\q\ is the sound mode, and b'{x) — 
— (3e^^ / {e^''^ — 1)^. The origin of this result is made 
clear by considering the analytic continuation of Eq. (pi) 
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to real frequencies, and summing the two terms. The 
pole at w = —Sq < (responsible for a finite contribu- 
tion at T = 0) is cancelled in favor of a double pole at 
Lu — Eq > 0, leading to Eq. (||), where the standard Bo- 
goljubov reduction of ps in a superfluid bosonic system is 
recognized. Thus, by fully including the dynamic struc- 
ture of the interaction for the phase, we obtain that in a 
Galilean- invariant system ps = p at T — |^ . 

A further consequence of the microscopic derivation 
of the effective action is a reduction of the strength 
x) of the static interaction term (V6')'* with respect 
to the XY model, as already observed in Ref. In 
the classical regime for PF only the term = of 
(i) and (H) should be considered. Therefore in Eq. 
(5) only the correction to D coming from the static 
interaction term oc (x/m^)(V0)^ survives, leading to 
6D Ki — Tx/d^m^-D^Q , qualitatively similar to the result 
of the classical XY model, SDxy ~ ~2T / (f S,^^"^ . Since, 
however, xjrn^D ~ we find that 



5D/5DxY ^ (kFCoY 
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i.e., in the classical limit SD is smaller than within the 
XY model, as far as kp^ < ^o- 

Let us now consider the effect of the Coulomb in- 
teraction between the electrons. At Gaussian level the 
density-density bubble xil) is dressed by the random- 
phase series of the Coulomb potential ^(q) = Ae^/|q|''^^ 
(here A is a constant depending on the dimension d and 
on the dielectric constant eoo), and the sound mode Sq 
of Eq. (^ is converted into the plasma mode uiq of 
the d-dimensional system. In deriving the anharmonic 
terms in Seff, we must now include the RPA density 
fluctuations in all the vertices. The one-loop correc- 
tions to D are formally identical to Eqs. (||)-(|1), with 
X — *■ Xlr- Thus at T = we recover again the can- 
cellations (^ of the bosonic diagrams, with Sq —> LUq. 
At the same time, since Xlr{<1,0) ~ l/^(q) vanishes 
as q ^ 0, the classical (cj„ — 0) term in Eq. (||) gives 
SDlb^ « -T/d{2d- l)Dm'^\e^^l'^~^, whereas the quan- 
tum XY model leads to the same result as the neutral 
case. Thus we estimate 



5Dlr/5Dxy ^ (Ef/Ec) {kp^o) 
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where Ep is the Fermi energy and Ec = Xkpe^ is a 
characteristic Coulomb energy. While within the XY 
model the Coulomb interaction modifies only the low- 
temperature behavior of ps , within the continuum model 
it affects also the high-temperature classical regime. 

To extend the previous results to a lattice model, we 
rewrite Eq. (^) introducing a hopping t between nearest- 
neighboring sites on a cubic lattice of spacing a, so that 
the free-electron dispersion is = —2t cos ak^ — p,, 

and we obtain the generalization of Eq. (p|) to the lattice 
case. Since we find that the cancellation of fermionic and 
mixed one- loop corrections to £) at T = still holds, we 



focus on bosonic corrections only. Differently from the 
continuum case, each insertion of a spatial derivative of 
9 in the fermionic loops is associated to a same-order k- 
derivative of ^k- Therefore, the (V6')^ term in Eq. (||) 
carries a factor = ^d^£,k/dk^. At Gaussian level, 
then, D{T - 0) ^ Ek,a ^^(1 - UEu), which, in the 
nearest-neighbor cubic model, equals the average kinetic 
energy, rather than p/m. In the anharmonic action of 
Fig la the factor A^ appears in each vertex with two in- 
coming 0-lines: as a consequence, the first two diagrams 
of Fig. lb, which give corrections to p in the contin- 
uum case, are now corrections to D, which are different 
from (and therefore do not cancel with) those coming 
from the shift of p, at fixed p. Moreover, on the lat- 
tice, the coupling of the fermions to 9 generates higher 
than second-order derivatives of 61 in (^) : in addition to 
the diagrams depicted in Fig. lb we must include also 
a diagram T, with four derivatives of 9 incoming in the 
same vertex, and a factor d'^£,]i/dk^, proportional to the 
average kinetic energy. The resulting bosonic one-loop 
correction to D is 

- E W i^XEEiq) + XEE 

q 

+ id-l)xEE + T -dx-'x%] 

- ujI {xIe^P^^I) - ^dVpEiq) + dx-'xpEV{q)] } ■ (9) 

The Xab bubbles correspond to the insertion of one (xpe) 
or two equal/different (xee/xee) factors A^; rjpE and rj 
are the bubbles with three Gq lines and one or no factor 
Aq respectively, rjpp corresponds to the first diagram of 
Fig. lb. The first terms in the two square brackets are 
the lattice analogous of the two terms of the continuum 
case, Eq. (^), whereas the last contributions come from 
the shift of p to keep p fixed. 

We perform the ujn sum in Eq. (^ with the PF Gaus- 
sian propagator P{q), and retain the leading order in 
C ^ l/$0 7 calculating the Xat bubbles at zero incoming 
momentum, and carefully evaluating the small q limit 
for rjpE and rj. We thus obtain an overall correction to D 
which is finite at T = in the lattice case. 



6D 



^^f^ [^Xee + {d- 1)xee + T ~2x, ^xIe 

I -1 2 1 9{d + 1) r -1 ^ / \ 

dx XpE] — [VpE - X XpEV\ ■ (10) 



At T = 0, X = nEk^k, with A^k - ^l/El; Xat = 

h Z^k^^ab^k, with TpE = Aa, TpE = A^ , tpE = AqA/3 

(a ^ /?); VpE = nEk^aMk with Mk = AUJE^, 
fj = ^Ek^k; ^ = a?D{T = 0). In the limit a 
0, t ^ oo, keeping 2a?'t = Xjm finite, Aq, tjm^ 
XEE, XEE -> x/f"-^ XpE xl^n, ijpE -> v/m, while 
T = a^D 0. Thus, Eq. (|l^) recovers the contimmm 
(Galilean-invariant) result 6D = 0. 
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We next turn again to the issue of the comparison be- 
tween SD and SDxy in the lattice case. At T = Eq. 
( p^ ) leads to an estimate of 5D/5Dxy of the same or- 
der of Eq. within a numerical factor. This is also 
true in the classical limit, both for the neutral and the 
charged system. As we discussed above, at high temper- 
ature only the static correction to D contributes, with 
a coefficient controlled in the neutral case by xee,Xee 
and T, which plays the same role (and has the same 
estimate) of x/"^^ i'^ the continuum case. The pres- 
ence of Coulomb forces does not change qualitatively 
this conclusion, and introduces only minor quantitative 
corrections, since the RPA expressions for xeEtXee 
have a finite limit for q — s- 0, contrary to the contin- 
uum case, e.g., x¥e = Xee - Xp£;^(q)/[1 + ^(q)x] - 
{xee-xIe/x) + xIeH'^^^ lx><e^- As a consequence, Eq. 

gives a proper estimate in the lattice case, even in the 
presence of long-range interactions. 

According to Eq. (|^) the PF-induced depletion of ps 
can be quite small at weak and intermediate coupling 
(particularly in the BCS limit) both in the quantum and 
in the classical regime, while it is of the order predicted 
by the XY model in the strong-coupling hmit kp^Q « 1. 
In this case the corrections to ps due to PF are size- 
able (even though less important than the quasiparticle 
contribution in determining the low-temperature depen- 
dence Q). The evaluation of Xab, fjpE for t/f/ <C 1 would 
lead to the conclusion that all contributions coming from 
the dynamic vertices are subleading with respect to those 
which arise from the static interaction X]Q=i(^a^)^ @- 
As a consequence, at strong coupling only the static in- 
teraction survives, analogously to what is assumed in the 
quantum XY model. Moreover, in the case of s-wave 
pairing, the value of the coefficient of X]a=i(^a^)^ ex- 
actly the same of the XY model (Q), with substituted 
by the lattice spacing a, in agreement with the strong- 
coupling expectation. For c?-wave pairing the coefficient 
of the static interaction changes only by a numerical fac- 
tor, in units of a, with respect to the XY model ]T^ . 

However the situation is more involved, specifically 
for the s-wave pairing symmetry. Let us consider the 
negative-?/ Hubbard model. In order to analyze the 
strong-coupling limit we need to include (i) the RPA fluc- 
tuations, also induced by U , in the particle-hole channel, 
and (ii) the fluctuations of |A| ~ [U /2)yJ p{2 — p), which 
fluctuates because p fluctuates |^ . When this analysis is 
carried out the contributions from the dynamical vertices 
are not subleading, and those from the static vertices do 
not reproduce the XY result by themselves. Neverthe- 
less, at T = the inclusion of both dynamic and static 
corrections to D leads again to SD — 6Dxy, with 
substituted by a. This result holds also in the classical 
regime, even though only approximately and provided 
the particle density is not too small [[l5[. 
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